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emphasis on the issues of reducibility and recoil corrections. An inconsistency when one pion ex- 
change is included in the production operator is demonstrated and then resolved via the introduc- 
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pv^ ' tion of "wave function corrections" which nearly vanish for static nucleon propagators. Inclusion 

of the recoil corrections to the nucleon propagators is found to change the magnitude and sign 
of the impulse production amplitude, worsening agreement with the experimental cross section by 
^ \ ~ 30%. A cutoff is used to account for the phenomenological nature of the external wave functions, 



and is found to have a significant impact for A < 2.5 GeV. 
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I. INTRODUCTION 

The theory of pion production in nucleon-nucleon colhsions has seen many advances. For 
an older review see Ref. l| and for a newer one see Ref. [2|. Modern calculations are 
based on effective field theory (EFT) in an effort to eventually obtain model-independent 
predictions. Such a goal is quite ambitious and naturally there are several major obstacles to 
overcome. Two such obstacles are the large threshold momentum and the presence of initial 
and final state interactions. Understanding pion production is interesting largely because 
of the aforementioned large threshold momentum. When neutrons and protons with such 
momenta collide, they come much nearer to each other than they do in the similar np — )■ d'j 
reaction. Thus we are able to "see" and therefore study a short-range region of strong 
nuclear force. 

The EFT that has been used to describe np — )■ dir^ is called baryon chiral perturbation 
theory |3|-l6|, a theory based on the low-energy symmetries of QCD. The theory consists of an 
infinite sum of interactions which are organized according to a "power counting" scheme, an 
expansion in the typical momentum appearing in an interaction divided by the symmetry- 
breaking scale A^ ^ ni]^. This power counting was developed assuming the momentum 
to be ~ ttItt- In order to produce a pion at rest, the initial relative momentum of the 
colliding nucleons must be p ~ ^Jm~m^ which means that the expansion parameter of the 
theory, x — v^"V"^/'"l/v = 0.38, is much larger than the normal 771^^/7711^ = 0.14. It was 
first proposed in Ref. \f\ that the terms in the theory be reorganized to reflect this large 
momentum. This modified power counting scheme is referred to as MCS. 

n 

In Ref. [8[, MCS was used to answer two important questions: why the theoretical 
pp — 7- (i7r+ (equivalent to 7ip — )■ dn^ under isospin symmetry) cross section was so much 
smaller than in experiment and why there existed formal inconsistencies with the next-to- 
leading-order (NLO) loops. The answers to both questions followed from a subtle reducibility 
issue with the rescattering diagram that dominates the cross section. In this paper we 
investigate the question of reducibility in the impulse approximation. 
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An attempt to address this issue was put forth in our recent study [9|] of charge symmetry 
breaking in 7ip — )■ rfvr^ where we introduced "wave function corrections." These corrections 
were calculated for the final state, and found to be a small fraction of the impulse diagram 
that they are correcting. However, this calculation suffers from a particular approximation 
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FIG. 1: Pion production operator. Solid lines represent nucleons, dashed lines represent pions, and 
ovals represent interactions. 

which we will describe. Fixing this approximation has a very significant effect. Furthermore, 
we show that the wave function corrections in the initial state are larger than one would 
expect in the MCS scheme. 

In Sec. [TTlwe review the np — > dn^ reaction and the impulse approximation's role. Then, 
Sec. mil examines the inconsistency that is found when one includes static one pion exchange 
(OPE) with the impulse approximation. Also in this section, wave function corrections are 
presented as a solution to the problem. Section [IV] discusses the correct implementation of 
the recoil corrections to the nucleon propagators. The effects of including a cutoff are shown 
in Sec. El Finally, we discuss the total cross section in Sec. |Vl]and conclude in Sec. IVIIi 

II. PION PRODUCTION 



The pion production operator in momentum space depicted in Fig. [T] is a function of the 
pion momentum q and I = k — p, where k (p) refers to the final (initial) relative momentum 
of the nucleons. The momentum transfer between the two nucleons in Ai is defined as 
q' = P2 — k2- We use a phenomenological, non-relativistic potential {V) that is static: the 
energy of each individual nucleon is conserved by V. Given this choice, working in the center 
of mass frame requires g' ° = ujq/2. It should be noted that this is an approximation (the 
"fixed kinematics approximation") that one needs to adopt in order to work in position space. 
If one works in momentum space and fixes q' ^ via energy conservation at the NNtt vertex, 
this is called the "equation of motion approximation." As was shown in Ref. lOj, both of 
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FIG. 2: Feynman rules 



these approximations have problems, particularly when considering initial state interactions. 
Nevertheless, it appears that the former is preferable to the latter. 

In this work we employ threshold kinematics, where q = (r?7,^,0), q' = {171.^/2,1), pi2 = 
(m7r/2, ±p), and ki^2 = (0,±A;) with \p\ = 359 MeV. Also, at threshold only s-wave pions 
(/tt = with respect to the deuteron) are produced and the initial state is purely ^Pi. We use 
the hybrid methodology introduced in Ref. 11| where "operators" (two-particle irreducible 
diagrams) are calculated perturbatively and then convolved with A^A^ wave functions which 
are obtained using phenomenological potentials. Details of how this procedure is carried 
out for A^A^ — )■ NNn can be found in Sec. Ill of Ref. J9|. Appendix [X] discusses the chiral 
Lagrangian that defines the theory. Using Eqs. ( 1A2|) and ( JA3I) . we obtain the Feynman 
rules shown in Fig. [2l 

At leading order, C(x^), the s-wave amplitude is dominated by the "rescattering" dia- 
gram, where a single pion is emitted from one nucleon and inelastically scattered by the 
other nucleon into an on-shell-produced pion. The only other leading order s-wave diagram 
is the impulse approximation, shown in Fig. [31(a) [3] . The details of the calculation of this 
diagram are given in Appendix [B For the initial and final states, we make use of three 
different potentials: Argonne vl8 12], Nijmegen II [l3|, and Reid '93 [l3|. The results for 
the s-wave reduced matrix element Aq are shown in Table [T] alongwith the leading order 
(LO) rescattering results which are detailed in Appendix C of Ref. |9(|. Aq also receives loop 
contributions at NLO (x^) which were calculated in Ref. [M]. However, a major result of the 
aformentioned Ref. [8i] was that these loops cancel, with the modification of the tree-level 
rescattering diagram being put on-shell. This modification (which is truly LO) is included 
in Table [B Finally, we point out that the experimental data (Sec. IVII) imply a reduced 
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FIG. 3: Impulse approximation operator alone (a), with OPE pulled out from the final state (b), 
and with OPE pulled out from the initial state (c). Solid lines represent nucleons, dashed lines 
represent pions, and ovals represent interactions. 

TABLE I: Reduced matrix elements of the rescattering and impulse production operators for three 
different potentials. 



Diagram 


Avl8 


Nijm II 


Reid '93 


/ires 
^0 


76.9 


83.4 


80.3 


.imp 
"^0 


4.9 


1.3 


3.5 



matrix element of 80 < Aq < 94. 

It is important to notice that the impulse approximation operator of Fig. [3]^a) cannot 
formally be convolved with the initial and final states as described above because the nucleon 
emitting the pion cannot remain on-shell. To put it another way, g"^ = because there is no 
way for the energy to be transferred. The common approximation made in pion production 
calculations is to ignore this formal difficulty. On the other hand, one "pulls" an OPE from 
the final state wave function [Fig. |3]^b)] in order to argue that the impulse approximation is 
leading order. This can be seen as an application of the Lippmann-Schwinger equation for 
the final state deuteron, 

I ^d) = GV I V'd), (1) 

where G represents the two-nucleon propagator and V represents the full potential. In Figs. 
^h) and[3]^c), we have made the approximation that V ~ OPE, which is known to be a 
good approximation for the deuteron 15|, but is not as valid for the initial state. Eq. ([1]) 



TABLE II: Reduced matrix elements of the impulse approximation with an OPE pulled out from 
the final state interaction [see Fig. [3]|^b)]. 



Diagram 


Avl8 


Nijm II 


Reid '93 


.OPE,red,f 
^0 


75.2 


64.6 


79.3 


.OPE,irr,f 
^0 


75.6 


64.7 


79.8 


-OPE,irr,f <OPE,red,f 
^0 ^0 


0.5 


0.1 


0.5 



begs the question: are the diagrams in Figs. M^a.) and Mjo) the same size? If so, we will be 
able to conclude that Fig. |3]^a) is doing its job well. 



III. INCLUDING OPE 



Calculation of Fig. Mjo) using Eq. ([T]) is detailed in Appendix O For lack of a better 
name, we will call this the "OPE reducible" diagram. In this calculation, we take G = 
{E — Hq)^^ = {— Efi—p'^ / 771]^)^^ where Ed = 2.22 MeV is the binding energy of the deuteron 
and 771^/2 is the reduced mass. The energy of the exchanged pion in this case is taken to 
be g'" = 0. This choice is consistent with the fact that the OPE is the first term in the V 
of Eq. ([1]) , which should be the same V that is used to generate the initial and final wave 
functions. The results are shown in the second row of Table Ull 

We find an inconsistency between the impulse approximation [Fig. |3]^a)] and OPE re- 
ducible [Fig. El^b)] diagrams: although they are equivalent according to the Lippmann- 
Schwinger equation, they are of very different size numerically. Using Avl8, they are 4.9 
and 75.2, respectively. Of course this inconsistency is not surprising when one notes that 
three-momentum transfer is provided for in the latter diagram but not the former. 

To resolve this problem, we re-consider the diagram in Fig. [3]^b) as a fully irreducible 
operator (without mentioning external wave functions). This is justifiable in that the left 
intermediate nucleon is off-shell by m,r, niore than the rri^/mAr typical of reducible diagrams. 
If we view it in this way, we are free to chose the energy of the exchanged pion to be 
g'° = m^/2 as mentioned in Sec. [TTl Additionally, the single-nucleon propagator for the left 
intermediate nucleon is taken from the rules shown in Fig. [2j The reduced matrix element 



for this "OPE irreducible" operator is given in Appendix [Cl and the results are shown in 
the third row of Table |lll We find that this diagram, which correctly accounts for energy 
transfer, is approximately equal to the OPE reducible diagram. 

The question remains: should Fig. El^b) be included, and if so, how? Until a clear proce- 
dure is defined for going from the full four-dimensional irNN coupled-channels formalism to 
the more common three-dimensional uncoupled formalism, this question is open to interpre- 
tation. We continue to take the view proposed in Ref. 9|, which is that the OPE irreducible 
diagram should be included with the OPE reducible diagram subtracted off to prevent dou- 
ble counting. This difference is shown in the fourth row of Table [ITl and is referred to as the 
(nearly vanishing) "wave function correction." Note that the cancellation is not as trivial 
as it appears. Schematically, OPE is two derivatives on a Yukawa that has different ranges 
for the reducible and irreducible cases. In going from reducible to irreducible, the radial 
integral gets bigger because the range increases. The derivatives bring down inverse powers 
of the range such that the overall amplitudes are similar in size. 

As in Sec. [TTl for s-wave pions at threshold we have a a ■ (a: + Pf) at the vertex where 
the pion is produced. For this reason the authors of Ref. 9] only considered OPE in 
the final state, assuming initial state OPE [Fig. [3t^c)] to be suppressed by the small final 
state momentum k. However, one needs to be careful when applying power counting to 
calculations that involve external A^A^ wave functions. At small distances the momenta of 
the nucleons (derivatives in position space) are distorted away from their constant values at 
asymptotically large distances. As an example of this difficulty, it can be shown that the 
k '^/2mN operator becomes larger than the p'^/2mN operator in the rescattering diagram. 
For this reason, we also calculate Fig. |3]^c) (for the details, see Appendix [C]). The results 
of this calculation are shown in Table IIIII where, again, the full wave function correction is 
nearly zero. 

There is one more formal point to discuss with regard to the above calculations. Expres- 
sions for nucleon propagators in irreducible diagrams differ based on the power counting 
scheme used. Consider the situation shown in Fig. HI Starting from the full, relativistic 



propagator, the authors of Ref. 16| showed that, if p° ~ ttItt and p ~ ^/rr^/rrT^, the correct 



propagator after emitting the pion is 

I (2) 

-gO + (2p • q- q"^)/2mN + ie ^ ' 



TABLE III: Reduced matrix elements of the impulse approximation with an OPE pulled out from 

the initial state interaction [see Fig. [3||c)]. 



Diagram 


Avl8 


Nijm II 


Reid '93 


.OPE,red,i 
^0 


-11.2 


-23.7 


-15.0 


/,OPE,irr,i 
^0 


-11.2 


-23.7 


-15.0 


/.OPE,irr,i <OPE,red,i 
^0 ^0 


~0 


~0 


~0 






P = miyv + p 

FIG. 4: Momenta of the nucleon propagator where v = (1,0,0,0) is used 

where the second term in the denominator is the recoil correction, suppressed by one power 
of X if 0'° ~ y/fT^TT^T^N- Note that we have ignored the resulting vertex corrections and 
antinucleon effects, both of which are suppressed for any choice of q^ . Equation ([2]), which 
we will refer to as the "new" method, comes in opposition to the "old" method which derives 
the propagator from the non-relativistic chiral Lagrangian, 

(3) 



pO _ qO _ (jj_ q^'^ jlvti^ -|- ie 

where the third term comes from the NLO Lagrangian and is a candidate for promotion in 
MCS counting. Let us now consider the case of Fig. El^b), where in terms of the momenta 
shown in Fig. Hlwe have p = {m.„/2,f>) and q = (m^,0). Since p^ — q^ ^ m.„, we promote 
the recoil corrections in the old propagator and find 



iG^:£ 



-m^ 



iG\ 



irr,f 
old 



m7r/2 — m^r — p"^ /2mN 
Next, consider the case of Fig. [3]^c) where we have p = {mT^/2,p} and q 
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(4) 



-ijj,p — k), and 



thus 

iG 



irr,i 
new 



u + (p2 — k'^)/2m]y + ie 



I 



^CoM = "-^ • (5) 

m-^jl + u — k '^ /2raN + ie 

Note that in the absence of distortions (|p| ~ y/n^/rriN, \k\ ~ 0), iG^ew = iGo\d for both the 
initial and final state propagators. For the sake of clarity we will define as the "free recoil 
approximation" (FRA) the use of these free particle values for the nucleon momenta. 

IV. NUCLEON PROPAGATOR RECOIL 

In Sec. Illlt the FRA was used for the recoil corrections to the nucleon propagators. In 
this section we calculate the diagrams again, treating the momenta properly as operators 
instead of numbers. In this section we use the old nucleon propagators for the irreducible 
diagrams. In doing so, we avoid the Gnew of Eq. ([5]), which would be difficult to evaluate 
exactly in position space. It was pointed out in Ref. [l7| that the old nucleon propagators 
have formal convergence problems owing to the large external momenta. Nevertheless, we 
expect to gain insight into the validity of the FRA using these propagators. 



iG 



irr,f 



-m^/2 — p 2/2m7v ' 



iG"'^' = ^ . (6) 

m^ — k'^/2mN + ie 

Of course, according to MCS, k should not be counted as ~ y/m-^niN and the k '^ /2m^ term 
should therefore not appear as in Eq. ([6]) until higher order. We choose to retain it here as 
an investigation into the effects of the distortions. For the reducible diagrams we continue 
to use G'""^ = {E - H^)-^ . 

The matrix elements can be calculated exactly in position space with Green function 
methods (see Appendix [D]) . The results are shown in Table HVl We find that the final state 
wave function correction evaluated without the FRA gives a sizable negative contribution 
of approximately —10. Additionally, we find that the initial state corrections become as 
important as the lower-order final state corrections. To verify the surprising results of this 
calculation, we examine as an example the ttim -^ C)0 limit of the radial integral for the 



TABLE IV: Reduced matrix elements for the wave function corrections with proper treatment of 
the momenta using the old expressions for the nucleon propagators. 





Avl8 


Nijm II 


Reid '93 


.OPE,irr,f 
^0 


80.8 


70.8 


89.1 


.OPE,red,f 
^0 


92.4 


81.2 


103.3 


.OPE,irr,f <OPE,red,f 
^0 ^0 


-11.7 


-10.4 


-14.2 


-OPE,irr,i 
^0 


5.1+24.5i 


15.1+34.7i 


7.8+27.8i 


4 0PE,red,i 
^0 


16.2+8.2i 


22.9+9.9i 


18.2+8.7i 


.OPE,irr,i ,OPE,red,i 
^0 ^0 


-ll.l+16.3i 


-7.7+24.8i 


-WA+19.li 
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FIG. 5: (Color online) Irreducible initial state OPE integral as a function of rrifq using Avl8. The 
solid line displays the FRA result and the real (imaginary) part of the exact propagator result is 
shown as a dashed (dotted) curve. 

irreducible initial state OPE in comparison with its analog from the previous section (which 
is independent of mj^i since /c = is used) , 

drr- (v^l^ + (I + ^) ^) G°™'^"'^(2/(mj2,r) + ,(mj2,r))i?.(r), 

(7) 

where the functions / and g are defined in Appendix O As shown in Fig. |5l in the large-mAr 
limit the recoil term in the propagator vanishes and we recover the leading order result. 
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V. CUTOFF DEPENDENCE 

It should not come as a surprise that odd things are happening in the short-distance part 
of the wave function, especially when we take derivatives. In the hybrid formalism we are 
using, this domain of the wave function is calculated from a phenomenological potential: 
Woods-Saxon for Avl8, one boson exchange for Nijm II, and Yukawa for Reid '93 (the very 
short range is exponential, exponential, and dipole, respectively). Because these potentials 
are fitted to experimental phase shifts, the wave functions derived from them can be consid- 
ered as infinitely high order in the EFT. Thus one should consider using a cutoff to account 
for this mismatch between the operator and the wave functions. Use of such a cutoff is 



referred to as EFT*, and was introduced in Ref. 11|. 

In this section we investigate the effects of cutting off the convolution integrals that 
account for the the presence of initial and final states as discussed in Sec. HVl We use the 



procedure of Ref. 18j, which modifies the Fourier transforms with a Gaussian cutoff, 

dH 



M{f^ = j^^/-Sl{l^)Mil) 



S. (/-) - exp J-ii j . (8) 

Note that the impulse approximation is not affected by such a cutoff scheme. For the OPE 
operators we define g\{u,r): 

lx{uj)gA{uj,r) _ r dH ^a.f._p,f,2 1 

47r J (27r)3 P + ii{ujf' ^ ' 

The exact evaluation of this integral and of the derivatives required to compute the diagrams 
of this work are shown in Appendix [El As desired, the cutoff regulates the behavior of g{r) 
at the origin, as shown in Fig. |6l The cutoff dependence of various reduced matrix elements 
is shown in Fig. [71 

The fact that we observe significant cutoff dependence of the wave function corrections 
above the typical scale of ~1 GeV is surprising. Indeed, this sensitivity indicates the need 

'or a counterterm because observables must be cutoff independent. As pointed out in Ref. 

19|, if one considers the difference of terms that comprise the wave function correction, 

1 1 ml/A 1 



g*'2-|-m| g*' 2 + 3772^/4 g/2_|_^^2 g'2 _)_ 3^72^/4' 

11 
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r(fm) 



FIG. 6: (Color online) Comparison of g (solid curve) and ^a (dashed curve) with lo = m-,^/2 and 
A = 1 GeV. 
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FIG. 7: (Color online) Cutoff dependence of various reduced matrix elements for Avl8 (solid curve), 
Nijmll (dashed curve), and Reid '93 (dotted curve). 

it can be argued that wave function correction is N^LO in the MCS scheme. However, if 
this view is to be accepted, the fact that the wave function corrections are much larger in 
magnitude than the LO impulse approximation should be considered surprising. 
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TABLE V: Reduced matrix elements for three different potentials. 



Diagram 


A^ 

A = oo 


18 

A = 1 GeV 


Nij 

A = oo 


mil 
A = 1 GeV 


Reic 

A = oo 


1 '93 
A = 1 GeV 


Rescattering (NLO) 


76.9 


76.2 


83.4 


81.5 


80.3 


79.1 


Impulse 


4.9 


4.9 


1.3 


1.3 


3.5 


3.5 


Final wfn cor (FRA) 


0.5 


-0.3 


0.1 


-0.7 


0.5 


-0.5 


Final wfn cor (exact) 


-11.7 


0.4 


-10.4 


-1.5 


-14.2 


-0.7 


Initial wfn cor (FRA) 


^0 


0.1 


^0 


0.1 


^0 


0.1 


Initial wfn cor (exact) 


-ll.l+16.3i 


-13.0+11. 7'i 


-7.7+24.8i 


-13.0-l-14.2i 


-10A+19.1i 


-13.0-hl3.0J 


Total (FRA) 


82.2 


80.8 


84.8 


82.2 


84.2 


82.2 


Total (exact) 


59.0+16.3i 


68.5+11.7i 


66.7-h24.8i 


68.4+14.2i 


59.2+19.1i 


69.0+13.0i 



VI. CROSS SECTION RESULTS 

Shown in Table |V] is a summary of the findings discussed in this paper at A = oo and 
A = 1 GeV along with the rescattering diagram. We also discuss the total cross section, 
which near threshold is parametrized as 



where q = mT^rj. At threshold, one can only calculate a, 



(11) 



m. 



a 



USn^sp' 



14 |2 

/in , 



(12) 



where s 



[m^ + md) 



Note that charged pion production is related to neutral pion 



production by isospin symmetry (breaking is expected to be small in the total cross section). 
This symmetry is the reason for the 1/2 present in the definition of a. The most recent 
experimental data are shown in Table IVI[ 

The theoretical total cross section as a function of the cutoff is shown in Fig. |H1 The theo- 
retical results include all diagrams up to Oix^). Thus theory can assign a rough uncertainty 
to the threshold cross section of 2 x ;y^ ;^ 30%. 
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TABLE VI: Experimental total cross section parameters 



Experiment 


a (/xb) 


/3(/ub) 


np -^ dvr" [20] 


184 ±5 


781 ± 79 


pp — > (i7r+ (Coulomb corrected) [21j 


208 ±5 


1220 ± 100 


pp — >■ (i7r+ (Coulomb corrected) [22] 


205 ±9 


791 ± 79 


Pionic deuterium decay [23] 


252+^1 


N/A 



200 f 
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(a) FRA nucleon propagators 
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0.5 1.0 1.5 2.0 2.5 

(b) Exact nucleon propagators 



A (GeV) 



FIG. 8: (Color online) Cutoff dependence of the total cross section. Avl8 (solid curve), Nijmll 
(dashed curve), and Reid '93 (dotted curve). 

VII. SUMMARY 

Before the findings of this work, the total theoretical cross section at A ^ 1 GeV [Fig. 
8(a)| was in agreement with the most recent experiment (fifth row of Table IVip at approx- 
imately the la level. Fixing the FRA approximation decreases the cross section, and if we 
stop here [Fig. 8(b)| , the agreement between theory and experiment becomes more tenuous 
(approximately the 2.5o" level). A second conclusion regards the MCS power counting, which 
dictates that |p| ~ ^Jrr^^m^ while \k\ ~ m^. We find that the wave function corrections 
of the initial state are of similar size to those of the final state once the FRA is removed, 
contradicting the previous sentence. There exists a contact term (an NNNNtt vertex) at 
N^LO along with tree-level diagrams proportional to the q low-energy constants (LECs) 
and two-pion exchange loops. Since all the LECs except the contact term are fixed by other 
data, it will be interesting to see if that contact term is of natural size. 
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We plan to investigate further into the concept of reducibihty. Specifically, we would like 
to define a clear procedure for deciding what to include in the impulse diagram. There has 
been a lack of consensus in the literature as to the inclusion of OPE, and if it is included, 
how that should be done. Understanding this issue is important, not only for calculation of 
the total cross section but also for p-wave pion production, since the leading contribution to 
j9-wave pion production comes from the impulse diagram. 
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Appendix A: Lagrange densities 

We define the index of a Lagrange density to be 

u = d+^-2, (Al) 

where d is the sum of the number of derivatives and powers of rriT,-, and / is the number of 
fermion fields. This represents the standard power counting for nuclear physics. The t/ = 
Lagrangian (with isovectors in bold font) is [7| 

£(0) = 1 (dnf - ^mlTT^ + NUdoN + ^N^ (r ■ c? ■ Vtt) iV + ..., (A2) 

where r and a are the pauli matrices acting on the isospin and spin of a single nucleon. The 
"+..." indicates that only the terms used in this calculation are shown. 

The u = 1 Lagrangian includes recoil corrections and other terms invariant under 
5^(2)^ X 5^7(2)^. 



£(1) = J—N^^^N 






+ ..., (A3) 



2m]\f 2mN 

where we use the values given in Table IVIIi Note that the terms with the Cj low-energy 
constants that appear at this order do not get promoted in MCS for these kinematics and are 
thus not used. Also, the terms with the di low-energy constants do not contribute to s-wave 
pion production. Finally, the NNNN contact terms, Cs,t, do not contribute because we 
are using a potential with a repulsive core [Ri{r)Rf{r) — )■ as r — )■ for /j = 1, // = 0]. 
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TABLE VII: Parameters used. 



m^ = 134.98 MeV 



gA = 1.32 MeV 



ruN = 938.92 MeV 



/^ = 92.4 MeV 



Appendix B: Impulse approximation details 



Evaluating the isospin matrix element 



(00 I ri,3 I 10) 



and using the vertex rule shown in Fig. [2], we obtain for Fig. [3](a) at threshold 



(Bl) 



00 I MLip,k) I 10 



9a rn^ 



2U2m 



ai-{p + k){2TTy6\p-k) 



N 



(B2) 



where Ai' = M-j ^2m^ 2m p 2m ^ = Ai/N . Since we are using position space np wavefunc- 
tions, we Fourier transform the matrix element with respect to I = k — p, which is identical 
to q' at threshold, 



dH 

{2ny 



e''-'{2'KY6-\p-k) = l 



(B3) 



Note that we group the p and k with their respective wave functions prior to performing the 
Fourier transform 

ai ■ {p+ k)^ai- {-i'^np + i^d)- (B4) 

Thus the full position space operator is 

(OO I M',{r) I lo) = -^H^^i ■ (^np - V,). (B5) 

To calculate the diagram with rescattering on the other nucleon, we consider how each part 
of the left side of Eq. ( 1B5I) transforms under 1^2. Since the strong part of the Lagrangian 
is invariant under isospin, M. is invariant. The initial isospin ket |1,0) is invariant as well, 

but |0, 0) ^ - |0, 0) and V -^ -V. Thus, 



00 I M',^^{f) I lo) = -2^|^^5^- (^n, - ^a 



(B6) 



LN 



The final spin-angle wave function is that of the deuteron, while the initial state for 
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s-wave pion production is solely ^Pi, 






(B7) 



where we have absorbed the unobservable (since there is only one initial channel available) 
phase into the definition of the matrix element. The spin- angle matrix elements are calcu- 
lated, 

dRAr) 



(/(f) ||5- ("^-^j \\i{r))=A7Ti 



Rfir) 



or j^ y ' ^ y ' Q^ 



, (B8) 



where Rf{r) = \/2u{r)/r + w{r)/r, Rf^2{r) = V^u{r)/r — 2w{r)/r, and Ri{r) = Ui^i{r)/pr. 
Using Eqs. f lB6P and fIBSp . we have the final result for the reduced matrix element, 



/limp 
^0 



K 



drr-mOO I ®(/(f) ||)-M(f)(||z(f))® | 10)) 



, .3 9A m^ ^ 
NSn — —- A, 



2/^ 2mN 
drr^ 



9^^ 2^ 

or r 



dRfjr) 
Or 



Ri{r) 



(B9) 
(BIO) 



Appendix C: Including OPE details 



1. Reducible OPE 



Taking just the q' terms at the OPE vertices. Fig. |3]^b) is given by 



M'{p, k) 



9A 



—I 



2f,j ^''^''"'''-"'w^^nm 



xri.3 






-Ed - p'^/ruN 



0-1 



2raN 



2p 



(CI) 



I — w^. Adding to this expression emission from the right nucleon and 



where /i(w)^ = m: 

approximating p^ = mT^rriN as discussed at the end of Sec. IIIIl we find 

(00 I M'{p, k) I 10) = ^;^ai ■ q'a, ■ q'— ^ ^ 



N 



(00 I M'{r) I 10) 



8/3 2m 
87r/3 2m^ 



q'^ + lj{Oy-Ed-m^ 
p{Of{Si2f{oo,r) + ai-a29{uj,r)) 



S -p, 



-Ed - m^ 



S ■ V(C2) 
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where the V acts on the initial np wave function, Su = 3cti ■ f(T2 ■ f — (?i ■ (T2 is the normal 
tensor operator, and 

-li{uj)r 



g{u,r) 



f{^,r) 



fi{u)r ' 



1 + 



-^(w)r 



(C3) 



come from the Fourier transform [see Eq. ( JB3I) ] of the pion propagator. Next, we evaluate 

(C4) 
where 



u[Lij,r) uir , 



wir) 



.g(cu,r) + 2V2^^^fico,r) 



w[u,r) Wir] 



r 
uir] 



igiu,r)-2fiu,r)) + 2V2^^fico,r). 



(C5) 



Thus, 



{f{r)\\iS^2fioo,r) + a,-a,g{uj,r))S-V\\t{f^) = 4m{^^^V2(^ + ^ 



w{uj,r) f d 1 



dr r 



Mr(^Q) 



and we finally arrive at the full reduced matrix element, 

g\ m^ /i(0)3 



^OPE.red.f 
^0 



-N 



L/(, 



00 



2/3 2mjv -Ed - m^ 
drr^ 



L/(0), 



u(u,r] r: I d 2\ w{uj,r) f d 1 
dr r I r \dr r 



Ri{r). (C7) 



2. Irreducible OPE 



Finally, as described in Sec. Illlt for the irreducible diagram we use (— m^/2 
p'^ /2mN)^^ ^ (—771^)"^ for the intermediate nucleon propagator and take u = mT^/2. 

.3 ^ ,,/^ /o\3 



4 0PE,irr,f 
^0 



2f^ 2mN -m^ 



(a 
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3. Initial state OPE 

For OPE in the initial state, the isospin matrix element is (00 | ti^sTi ■ T2 | 10) = 1, and 
because the initial state consists of just one channel, ^Pi, 

(Sufico, r) + <?! ■ <T2 g{cu, r)) 1^ P,) = {2f{cu, r) + g{co, r)) 1^ P,). 

Evaluating the S ■ V reduced matrix elements, we find 

/,OPE,rcd,i _ AT ^A ^'^ ^^^^ /"^ Ti(f\\ 

.OPE,irr,i _ T^j 9a ^tt At(^7r/2) /"J jif^ ,^n 



{2f{u,T) + g{u,T))R,{T). (C9) 



U{u) = / rfrr^ V2- ^ + + _ ^ 

J \ or r \or rj r J 

Appendix D: Exact ^vave function corrections details 

Consider the nucleon propagator for reducible OPE in the initial state. Pulling out a 
-rriTv and expanding this function in spherical coordinates, we have 



-2im 



gig|r-r'| 



47r|r — r'\ 
= 2m^eE^''(^^<)^!'^(^^>)^- (^')>l(r), (Dl) 

l,m 

where ^ = ^JnQmjq and r<(r>) is the lesser (greater) of |r|, \r'\. This spherical partial wave 
expansion was derived from the differential equation 

where iG^ = —imj^Q. First, one solves the homogenous equation and requires both finiteness 
at the origin and outgoing wave behavior for large r. Thus, G{r,r') = Aj/(^r<)/iJ \^r^). 
Next, the boundary condition at r = r' is obtained by integrating the differential equation 
across the boundary. In terms of g{r,r') = rr'Q{r,r'), 

■^9>{r,r%=r'+e - ■^g<{r,r%=r'-e = -1, (D3) 

or or 



19 



which yields A = i^. At this point in the diagram, the two-nucleon state is still ^Pi, so we 
preform one of the angular integrals and obtain 



Wi 



Thus, 



^: 



OPE,red,i, 
' 



m-n/2) 



iGoir,r') -^ mNUi{Cr^)h\ ^(^r>). 



^^ m^MO)! (_,^^^^.(o)) 



(D4) 



L\mj2) 



VI 2m^ 3 
dr dr' r'^r''^ I \/2 



d u{r) 
dr r 



d 3\ w{r] 
dr r ) r 



xji(er<)/iS'Her>) (2/(m^/2,r') +^(m^/2,r'))i?.(r('p5) 



For the irreducible initial state OPE, the only difference is that a —2mN gets pulled out 



and the momentum becomes ^' = ^/2rri~rri^, 



A 



OPE,irr,i, 
> 



mj2) 



-N- 



{-2imNi'L\mj2)) 



(D6) 



2/3 2m^ 3 

For the final state OPE, we can obtain the correct Green function from Eq. ( ]D4|) by letting 
^ — 7- z^ and using the correct / for the term under consideration. 



Appendix E: Cutoff details 



In this section we display the exact expressions needed to implement the Gaussian cutoff 
of Sec lYl For the OPE diagrams, the integral of Eq. ([9]) is evaluated. 



9A{^,r) 



2 



^-- — erfc h 

ji{u)r \ 2 A 

^^erfc h ^-V^ 



(El) 
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One also needs derivatives of Eq. (lEl 






(SufAiuj, r) + ai- a2/A(w, r)) 



r 



M(a;)VA2 



A 



+ (yU — > yU and A -^ —A) 



erfc (-— + ^ 

fi{u)r {fi{u)ryj V 2 A 

A2 3 

yu(w)r + IH ^^ I e 



fi{uj)r 



Ar I mMN^ 
" 2 "•" A J 



/A(c^,r) 



2^' 



M(a;)VA' 



erfc(-- + ^ 



A 



A2 



+ {fi -^ fi and A -)■ -A) 



/i(a;)r + 1 e '^ 2 + a J 









(E2) 
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